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In this work we study the constraints on the dark matter interaction with the standard model 
particles, from the observations of dark matter relic density, the direct detection experiments of 
CDMS and XENON, and the indirect detection of the p/p ratio by PAMELA. A model independent 
way is adopted in the study by constructing the effective interaction operators between dark matter 
. . . and standard model particles. The most general 4-fermion operators are investigated. We find that 

T-H I the constraints from different observations are complementary with each other. Especially the spin 

. independent scattering gives very strong constraints for corresponding operators. In some cases 

' the indirect detection of p/p data can actually be more sensitive than the direct detection or relic 

CN ' density for light dark matter (< 70 GeV). 

Q ■ PACS numbers: 95.35.+d, 95.30.Cq, 95.85.Ry 
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I. INTRODUCTION 

The existence of a significant component of nonbaryonic dark matter (DM) in the Universe has been well confirmed 
^ • by astrophysical observations [ll-Q in recent years, however the nature of this substance remains unclear. Since there 
^ ' is no candidate for DM in the Standard Model (SM) of particle physics, it implies the existence of new physics beyond 
. the SM. Probe of the microscopic identity and properties of DM has become one of the key problems in particle 

physics and cosmology (for reviews of DM, sec, for instance, 
, Among a large amount of theoretical models, a well-motivated candidate for DM is the weakly interacting massive 
J> ' particle (WIMP). This WIMP must be stable, nonrelativistic, electrically neutral, and colorless. If the mass of WIMP 
CsJ \ is from a few GeV to TeV and the interaction strength is of the weak scale, they can naturally yield the observed relic 
■ density of DM, which is often referred to as the WIMP miracle 0] . A huge variety of new physics models trying to 
] solve theproblems of the SM at the weak scale can naturally contain WIMP candidates, such as the supersymmetric 
! models [l[i-[ni, extra dimensional models [H-Iiil, little Higgs models [lO-Sl, left-right symmetric models [UHlil, 
04 ' and many other theoretical scenarios. 

The above mentioned specific models are well-motivated, however they still lack experimental support. We do not 
know whether nature really behaves like one of them or some other yet unconsidered theories. Moreover, in case the 
DM particle is the only new particle within the reach of LHC and other new particle species are much heavier than 
DM, it will be very difficult to tell which model the DM particle belongs to. Additionally, it is possible that the 
DM may be first observed by direct or indirect detection experiments. These early observations may only provide 
information about some general properties of the DM particle, and may not be able to distinguish the underlying 
' theories. Therefore, the model-independent studies of the DM phenomenology arc particularly important for they 
may avoid theoretical bias [29l - l3^ . Recently there have been quite a few papers follo wing su ch consideration and 
adopting a model-independent way to study various phenomenologies related with DM p35l - [43j . Especially the relic 
density measured by WMAP direct detection from CDMS [ij], XENON [4^ and possible collider signals from 
LHC are considered in these studies. 

In this work, we first construct the general effective 4-fermion interaction operators between DM particles and the 
SM particles, which extend the effective fermionic WIMP interactions given in Ref. [s^l- Here we focus on Dirac 
fermionic DM. Discussions on scalar and vector DM will be presented in companion papers. We then give updated 
constraints from the DM relic density within the 7-year WMAP data @ and the spin-independent WIMP-nucleus 
elastic scattering searches by CDMS II [31 and XENONIOO [i^, and compare our results with those in ReL [3^ . 
In addition, we present new phcnomcnological constraints on these effective models from the spin-dependent WIMP- 
nucleus elastic scattering searches by CDMS [i^ and XENON i47j] and the cosmic-ray antiproton-to-proton ratio by 
PAMELA [sl]. We find that the constraints from different kinds of experiments are rather comparable. Combination 
of these constraints provides more information of the effective models. 
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This paper is organized as follows. In Sec. [IT] we briefly describe the effective DM models of various 4-fermion 
interaction operators. In Sec. Illli IIVI and |V] we explore the constraints on these models from the DM relic density, 
direct and indirect detection searches, respectively. In Sec. IVII wc discuss the validity region of effective theory and 
present the combined constraints on the effective coupling constants of these models. Sec. IVIII is the conclusion. 



II. EFFECTIVE MODELS 



We start with the case that DM is a single Dirac fermionic WIMP (x). Instead of considering a WIMP candidate 
from a specific theoretical model, we study the phenomenologies in a model-independent way by constructing effective 
interaction operators between x ^-nd the SM particles. These interaction operators are constrained only by the 
requirements of Hermiticity, Lorentz invariancc and CPT invariance. 

To proceed, we make the following assumptions similar to those in Ref. (s^: (1) The WIMP is the only new 
particle at the electroweak scale, and any new particle species other than the WIMP has a mass much larger than the 
WIMP. This implies that the WIMP's thermal relic density is not affected by resonances or coannihilations, and this 
makes it possible to describe the interaction between the WIMPs and the standard model particles in terms of an 
effective field theory. (2) The WIMP only interacts with the standard model fermions through a 4-fermion effective 
interaction, but not with other particles like gauge or Higgs bosons. This 4-fcrmion effective interaction is assumed 
to be dominated by only one form (scalar, vector, etc.) of the set of 4-fermion operators for simplicity. (3) The 
WIMP annihilation channels to the standard model fermion-antifermion pairs dominate over other possible channels. 
In other words, the possible channels to final states that include gauge or Higgs bosons are assumed to be negligible 
for simplicity. However, it should be noted that the supersymmetric DM model actually cannot satisfy the three 
assumptions above simultaneously in order to give the correct relic density. Even so these assumptions are still useful 
for a general research. 

The effective Lagrangian between two fermionic WIMPs (x and x) two standard model fermions (/ and /) is 
given by only one of the following expressions: 
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— 7=^X7''(l + 75)x/7/i(l + 75)/ (12) 
/ ^2 

Left handed-right handed interaction (LR) : Clr = ^ — ^^^X7^(l ~ 75)x/7ai(1 + 75)/ (13) 

G 

Right handcd-lcft handed interaction (RL) : Crl = Y^ — ^X7^(1 + 75)x/7a'(1 - 75)/ (14) 

where the sum of / is over all the standard model fermions, and the effective coupling constants G are real numbers 
which have mass dimension of —2. The 4 chiral interaction operators, £lLi 'Crr, £lr f^nd £rl, are just the combina- 
tions of the other operators mentioned above. Here we do not include the interaction operators involving derivative 
5^ insertions into fermion bilinears, for they have higher momentum dimensions and may be safely ignored in the 
small momentum limit. Note that the alternative tensor interaction term in (|10p has two other equivalent forms, i.e., 

Each form of the interaction operators listed above represents an effective model of the WIMP coupled to the 
standard model fermions. For each case, we can calculate the corresponding annihilation and scattering cross sections, 
which depend on the WIMP mass and the coupling constants G/. Associated with the recent results of the DM 
relic density, direct and indirect detection experiments, we can obtain the phenomenological constraints on G/. It 
would be interesting and meaningful to compare the constraints derived from different experiments. 

It is worthwhile to note the symmetry properties of these operators under discrete C, P, and T transformations. 
The first 5 forms of the operators, £sj £v, and £t are separately invariant under C, P, and T, while the 
transformation properties of the other operators are summarized in Table ID The transformation properties of the 
operators under CP are the same as those under T, given that the coupling constants G are real-valued numbers. 
Thus, all the operators are actually CPT invariant. If the future experiments indicate that there were some of the C, 
P, and T symmetries in the DM sector, we may use this table to concentrate on or exclude some interaction operators. 



TABLE I. The transformation properties of the 4-fermion operators under C, P, and T. Since £s, £p, £v, £a and £t are 
separately invariant under C, P, and T, they are not listed below. The plus means being invariant under the transformation, 
while the minus ' — ' means sign reversal. The transformation properties of the operators under CP are the same as those under 
T, given that the coupling constants G are real numbers. 
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III. WIMP ANNIHILATION AND RELIC DENSITY 



In order to determine the relic density of WIMPs and the source function of cosmic-ray particles produced by the 
DM annihilation in the Galaxy, which is relevant to the DM indirect detection, we need to calculate the cross sections 
of WIMP-antiWIMP annihilation to fermion-antifermion pairs for each case listed in the last section. The result is 
given by 
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where s is the Mandelstam variable, M^ is the WIMP mass, the sum is over the final state fermion species /, and 
Cf are the color factors, equal to 3 for quarks and 1 for leptons. The annihilation cross sections of the four chiral 
interactions, aLL, ann, cr_R_R, ann, CTL/j^ann and auL.ann, havc the exactly same formula, so they can be denoted by a 
common symbol aim 

with the corresponding coupling constants denoted by Gcj- Note that E cis.(IT5l) -([T7 )) agree 
exactly with Eqs.(6j~(8) in Ref. [HI, while Eqs.(ITHl),(IIll) are slightly different from' Eqs.(9),(10) in 

In the very early Universe, the WIMPs were in thermal equilibrium. As the Universe expands, the WIMPs departed 
from thermal equilibrium when they were nonrelativistic, and finally froze out to yield a cold relic roughly when the 
annihilation rate dropped below the Hubble rate. This evolution process is described by the Boltzmann equation 
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where H = d/a ^ ^87rp/(3M|[) is the Hubble rate with Mpi denoting the Planck mass, (n^) is the number density 
of WIMPs (antiWIMPs), and the thermal average ((7aiin'yM0i) will be explained below. For Dirac fermions without 
particle-antiparticle asymmetry, n,.^ = n^, and the total DM paritcle number density is tidm = 2n^ [i,|4i,|4|. At 
the early time, when the temperature T 3> M^, the WIMP number density was very close to its equilibrium value 
oc T^, and the annihilation rate per unit volume F = n^n^ (cTannVM^i) was much greater than the Hubble expansion 
rate per unit volume 3Hn^ and enough to maintain the thermal equilibrium. However, as the temperature T decreased 
below M^, the equilibrium number density was exponentially suppressed, n'^ ~ g[M^T/{2TT)]^/^ exp {—M^/T), where 
p = 2 is the number of degrees of freedom of a fermionic WIMP. Eventually, the annihilation rate became smaller 
than the expansion rate, and the WIMPs froze out of equilibrium. 

The thermally averaged quantity (cTann^'Moi) should be treated carefully. As pointed out by Ref. [i^, the M0ller 
velocity umqi in Eq. (|26p is defined by wmhi = \/ {pi ■ P2)^ — "z?™!/ (^1^2) = y/Wi ^ — |vi x V2P with subscripts 
1 and 2 labeling the two initial DM particles and particle velocities = Pi/Ei [i — 1,2). The M0ller velocity VM<t\ 
equals the relative velocity Wrei = |vi — V2I only when the collision is collinear, vi x V2 = 0. Since the Boltzmann 
equation, Eq. \2Q\ . is expressed in the cosmic comoving frame ji^, in which the gas is at rest as a whole, the thermal 
average {(t aYinVyi0\) must be taken in this frame. Fortunately, even including relativistic effects, it has been shown (ioj 

that ((TannWM0i) = (o'ann'yiab) '^'^ , whcrc Wlab = |vi,iab — V2,iab| and the right-haud side is computed in the lab frame, 
in which one of the two initial particles is at rest. Thus, it is convenient to calculate the thermal average in the lab 
frame using the method described in [4^ . 

Cold DM requires that the freeze-out of WIMPs occurred when they were nonrelativistic. In the nonrelativistic 



5 



limit, we can parameterize a, 
[isl |49| , we then obtain ((Tann 



annw = + 6^2 + 0[v^) 0, where u = viab = y^s(s - 4A/2)/(s - 2M2). According to 
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b in the effective models. Due to the common factor (s — 4M^) in Eqs. ([TS|) - (^5]) . s must be expanded up to 
order v'^ to get the correct coefficients h. In the lab frame, s = 2M2(1 + 1 /Vl - w^) = 41/^ + Af^w2 _^ 3^2^4 _^ ^^^e^^ 
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25|) and expanding (Tanni' in powers of v up to order u^, we obtain 
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from which, one can easily read off the corresponding thermally averaged quantities (canni')- Our results of 
Eqs. jSO]), dSH), and dST]) agree well with Eqs. (35) and (39) of Ref. Note that Eq.^ is the same as 

Eq.(13) in Ref. while the 0{v^) terms of Eqs.dMJ-lIM: and both the 0{v^) and 0{v'^) terms of Eq.^ are 

different from the corresponding terms of Eq.(14)-(17) in jS^]. In spite of these differences, they do not have much 
effect on the main results of Ref. [s^] . This is because the calculated relic density depends mainly on the leading 
term of (o-ann^^)- as we will see below. 

Using the standard procedure [1,01 to approximately solve the Boltzmann equation ([26]), we obtain a relic density 
of DM particles as 
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^ Note that this expansion in powers of v'^^^ is equivalent to that in Ref. [4^: (Tann''iab = o-'"' + a'-'^'e + C(e^) with e = (s — AM^) / {4Ml) . 
Since i^i^b = 2-^e(l + e)/(l + 2e), one easily obtains the relation between these two expansions: a'''' = a, a'^' = 4b, etc. 
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where Xf = My^/Tf with Tf being the freeze-out temperature, g*{Tf) is the total number of effectively relativistic 
degrees of freedom at freeze-out, Tq = 2.725±0.002 K [501 is the present CMB temperature. The freeze-out temperature 
parameter x j can be evaluated by numerically solving the following equation: 



where c is an order one parameter defined by the freeze-out criterion and determined by matching the late-time and 
early-time solutions. The precise value of c is not so significant for the numerical solution of Xf due to the logarithmic 
dependence in Eq. (P^ . and we take the usual value c = 1/2 in the calculation. Noting that in Eqs. and ([55)1 
depends on the temperature T, we adopt the recent numerical result of g*{T) in Ref. (5l| where the quark-hadron 
transition temperature is taken to be 200 MeV. 

According to the observed DM relic density, r^DM^-^ = 0.1109 ±0.0056 we can estimate the relation between the 
effective coupling constants G/ and the WIMP mass in each effective model of 4-fermion interaction operators, 
as shown in Fig. [T] Two kinds of coupling constants are considered here. In the upper frame of Fig. [TJ we show the 
results for the case when the effective couplings to all the standard model fermions are equal (universal couplings). 
In the lower frame of Fig. [1] we show the results for the case when the coupling constants are proportional to the 
fcrmion mass to/. In both cases, G/ decreases as increases for fixed Q.-oyih'^ in each effective model. Besides, 
Fig. [l]has several interesting features: 

• In both cases, the 4 curves of G/ vs. for the effective models of scalar, scalar-pseudoscalar, axialvector and 
axialvector-vector interactions lie well above the other curves. This comes from the fact that (Ts.sLim^^ ^sp.aun^ 
and (JAV. annV are of order 0{v'^)\ although the leading term of a a, a,nnV is of order 0{v^), it is smaller by a factor 
of mj/M-^ than the 0{ifi) terms for other types of interactions. 

• In the case of G/ cx nif, there is an obvious downward bend in the curve of G/ vs. at about ^ 
nit = 171.2 GeV in each effective model. This can be easily explained as follows. In the low velocity limit, 
the threshold for the annihilation channel x X ~^ ^ ^ is about ^ rrit- Since G/ (x m/, the WIMP couples 
much more strongly to the top quark than to other fermions, and the corresponding channel X X ~^ ^ ^ gives a 
tremendous contribution to the total (a'ann''^)- This finally makes the curve to bend down. 

• In the case of universal couplings, there are 5 pairs of nearly identical curves, because in each pair their 
corresponding fTannW differ only by terms of 0{v'^) and/or terms of m'j/M'^. This feature is also noted in 
Ref. [1^. These pairs are: (1) the curves for scalar and scalar-pseudoscalar interactions, and we denote this 
approximate identity of the two curves by S ~ SP for short here and henceforth; (2) P ~ PS; (3) V ~ VA; (4) 
T ~ C; (5) A ~ AV except for some small regions. From these approximate identities for pairs of curves, we 
see that the predicted relic density relies mainly on the leading term of (aann'^) in each model, especially on the 
C(M^) term in the leading term. 

• In the case of G/ cx m/, there are also 5 pairs of nearly identical curves except for some small regions. These 
pairs are the same as those in the above case, though deviations of the two nearly identical curves in each pair 
become large in this case. 

In Fig. [2] we show the curves of JIdm'i^ vs. for fixed coupling constants in the models of scalar, pseudoscalar, 
vector, axialvector, tensor and alternative tensor interaction operators. Due to the nearly identities described above, 
we do not include the curves for scalar-pseudoscalar, pseudoscalar-scalar, vector-axialvector, axialvector-vector and 
chiral interactions. Here we still consider the two kinds of coupling constants. In the upper left and upper right frames, 
we show the results when the coupling constants are proportional to the standard model fermion masses, G/ cx to/, 
and various values of the couplings are taken: G/ x (1 GbV/to/) = 10"®, 10"^, IQ-^, 10"^ and 10"* GeV"^ 
This proportionality of the couplings to the fermion masses may come from Yukawa couplings of a Higgs mediated 
interaction or some other unknown underlying mechanism. In the remaining 6 frames, we show the results for the case 
when the effective couplings to all the standard model fermions are equal (universal couplings), and various values of 
the couplings are taken: G/ ~ 10~^, 10^^, 10~®, 10"^ and 10^"* GcV^^. The curves in Fig. [2]bend down to more or 
less at about ^ 1.72 GeV, 4.20 GeV and 171.2 GeV, which exactly correspond to the masses of charm, bottom, 
top quarks, respectively. 

Comparing with the results in Ref. [32| , we observe that the curves in the first 6 frames of Fig. [5] are a little 
higher than those given by Ref. [32|. This slight difference may be caused by the following reasons: (1) We use 
JIdm = + = 2f2^ for Dirac fermionic WIMPs with the assumption of no particle-antiparticle asymmetry. (2) 
Some formulas of CannW we obtained differ slightly from those given by (3^ . as already described at the bottom of 
Eg. ([571) . (3) We use the effective degrees of freedom g*{T) given by |51| . 





(39) 




7 



10" 



10" 



— 1 1 1 1 1 — I — I — I — [— 



c<r 10"° 

I 

> 

O 

c3" 10"^ 



10" 




10" 



10' 



T — I — I — I 1 1 1 1 1 r- 

Scalar 

Pseudoscalar 

Vector 

Axialvector 

Tensor — ■ — 
Scalar-Pseudoscalar - - - - 
Pseudoscalar-Scalar - - - - 
Vector-Axialvector - — — 

Axialvector-Vector 

Alternative Tensor — ■ — 
Chiral 



Dirac Fermionic WIMP 
Fixed by Relic Density 



10^ 



(GeV) 



10^ 



10^ 



10" 



10" 



(M 

> 10" 

o 

10" 

E 



1 10" 



-8 



10 



10 



10 



11 



10' 



Scalar 
Pseudoscalar 
Vector 
Axialvector 
Tensor 
Scalar-Pseudoscalar 
Pseudoscalar-Scalar 
Vector-Axialvector 
Axialvector-Vector 
Alternative Tensor 
Chiral 




Dirac Fermionic WIMP 
Fixed by Relic Density 
Gf rPf 



—I I I I I I i_ 



—I I I I I I i_ 



10^ 



My (GeV) 



10' 



10^ 



FIG. 1. The predicted coupling constant Gf as a function of the WIMP mass M^, fixed by the observed relic density, 
^^dm'i'^ = 0.1109 ± 0.0056 in each effective model of 4-fermion interaction operators. In the upper frame, results are given 
for the case when the effective couplings to all the standard model fermions are equal (universal couplings). In the lower frame, 
results are shown for the case when the coupling constants are proportional to the fermion mass m/. In both cases, several 
pairs of curves are nearly identical. See the text for more details. 
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FIG. 2. The predicted thermal relic density (dashed lines) of Dirac fermionic WIMPs with scalar, pseudoscalar, vector, 
axialvector, tensor and alternative tensor interactions respectively. In the upper left and upper right frames, results are given 
for the case when the coupling constants are proportional to the fermion mass nif, Gf x (1 GeV/m/) = 10~*, 10~^, 10~^, 
10~^ and 10-* GeV-^ In the remaining 6 frames, results are shown for the case when the couplings to all the standard model 
fermions are equal (universal couplings), Gf — 10"*, 10"^, 10"®, 10~^ and 10~* Ge-V"^. The horizontal solid band shows the 
range of the observed DM relic density, fioivi/i^ = 0.1109 ± 0.0056, measured by WMAP 0]. 
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It is important to note that the results in Figs. [T] and [2] are found under the assumptions presented in Sec. |TT1 If 
resonances, coannihilations, or annihilations to final states other than fermion-antifermion pairs are significant, the 
actual curves in Figs. [T] and [2] will be significantly lower than these shown there, as pointed out in (32j |. 



IV. DIRECT DETECTION 



In this section we discuss the direct detection constraints on the effective models of Eqs.(fT|)- ([T4| . Direct detection 
experiments are designed to measure the recoil energy of the atomic nuclei when the WIMPs elastically scatter off 
them. The WIMP-quark interactions in the effective models naturally induce the WIMP-nucleon interactions, and 
the latter further induce the WIMP-nucleus interactions. Such interactions may lead to the elastic scattering of the 
WIMPs with the nuclei, which may be detected at the direct detection experiments. 

The velocity of the WIMP near the Earth is thought to be of the same order as the orbital velocity of the Sun, 
ti ~ 0.001c. Because of this small velocity, the momentum transfer in the WIMP-nucleus scattering is considerably 
small compared to the masses of the WIMP and the nuclei. Thus all the WIMP-nucleus cross sections can be 
calculated in the limit of zero momentum tranfer. In this limit the WIMP-quark interaction operators, £p, £sp, 'Cps, 
>CvA, 'Cav a-iid £rp, and their correspondingly induced WIMP-nucleon interaction operators have no contribution to 
the WIMP-nucleus cross sections and thus they are not sensitive to direct detection experiments. This is because 
in the zero momentum transfer limit some fermion bilinear operators become zero, for example, the operator tpj^ip 
vanishes, and the time component of ipj^jsip ^^nd the space components of ipj^ip vanish as well. For more details on 
this issue, see Ref. [s^ . 

Among the remaining operators relevant to direct detection, the scalar and vector interaction operators, Cs and 
Cy, are referred to as spin- independent (SI) interactions, while the axialvector and tensor interaction, £a and Ct, 
belong to spin-dependent (SD) interactions. And the chiral interaction operators are the combinations of SI and 
SD interactions. The SI interactions of all the nucleons add coherently in the target nucleus, and the corresponding 
WIMP-nucleus cross section is proportional to the square of the atomic mass number of the nucleus. On the other 
hand, since the spins of nucleons in a nucleus tend to cancel in pairs, the SD interactions rely mainly on the spin 
content of one unpaired nucleon and the corresponding cross section is not enhanced for heavy nuclei. 

We would like to illustrate the calculation in the effective model of scalar interaction operators. Eq.(IT]) can induce 
the effective Lagrangian for the WIMP-nucleon couplings, which reads 



E 

N—p^n 



Gs. 



N 



V2 



XXNN 



(40) 



where the WIMP's effective Fermi couplings to the nucleons (protons and neutrons), Gs,n {N = p, n), are related to 
the coupling constants to quarks by 



Gs,N — ^ Gs,qfq 



n'^n 



q—u^d.s 



E Gs,qfQ 



* q=c,b,t 



m„ 



(41) 



where the nucleon form factors are fP = 0.020 ± 0.004, = 0.026 ± 0.005, ff = 0.118 ± 0.062, = 0.014 ± 0.003, 
/? = 0.036 ±0.008, = 0.118 ± 0.062 [H, and 



fN _ 

JQ - 27 



(42) 



From Ea. (|40)) and its induced WIMP-nucleus interactions, it follows that the cross section for a WIMP (x) scattering 
elastically from a nucleus (A) in the zero momentum transfer limit is given by [s^ 



Scalar Int. 



TT (Mx + ttia)^ 



(43) 



where uia is the target nucleus mass. Z and {A — Z) are the numbers of protons and neutrons in the nucleus. The 
factor 1/2 in the square bracket comes from the fact that a Dirac fermionic WIMP and its antiparticle are different 
[5^ . However, this factor is missing in Eq.(20) of Ref.fs^]. Indeed, such a factor 1/2 does not exist in the expression 
for a self-conjugated WIMP such as a Majorana fermion, but it seems that a Dirac fermionic WIMP is considered in 
Ref. [3^ ■ otherwise the fermion bilinears X7^X and X'^'^'^x vanish for a Majorana fermion x- Taking the special case 
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when the nucleus is just the nucleon (proton or neutron) in Eq. (|43)) . we obtain the WIMP-nucleon cross section in 
the zero momentum transfer Umit: 



Scalar Int. 



■ mN) 



Gs.N 

~7T 



(44) 



Likewise, we compute the WIMP-nucleus cross section in the effective models of vector, axialvector, tensor and 
chiral interaction operators, resulting in 



Vector Int. 
Axialvector Int. 
Tensor Int. 
Chiral Int. 



ctt, xA ■ 

OC,xA 



Mlm\ 



7r(Mx + rriAY 



'x"'A 



Ja + 1 



16M^m\ J A + 1 
7r(Mx + mAY Ja 



GA,n „A 



Ga,p nA 

Gt.p ^a , Gt,,i „a 



V2 V2 



(45) 
(46) 
(47) 



Z 



.G 



V,p 



V2 
Ja + 1 



{A-Z) 



Gv,n 



1 2 



V2 



it{M^ + ruAf J A 



Ga,p gA _|_ GA,n ^^A 



^/2 " V2 



TT{M-^+mAY 



Z^ + ^A-zf""- 



V2 



V2 



(when A > 1) 



(48) 



where J a is the nuclear spin, is the expectation value of the total spin of the nucleon (iV) in the nucleus {A), and 
the WIMP's effective Fermi couplings. Gat, to the nucleons are related to those to quarks, Gg, by 



Gv.p — 2Gv^u + Gv,d , 
Gam = ^ GA,q^q 

q—u.d,s 

Gt,N — ^ GT,q^q 
q—u.d,s 



G\ 



Gv,u + 2G 



N 



(49) 
(50) 

(51) 



with the form factors = 0.842 ± 0.012, = -0.427 ± 0.013, Af = -0.085 ± 0.018 K = ^d^ ^d = ^u, 
A" = Af. In Eq. ipS]) the WIMP-nucleus cross section for chiral interactions receives contributions from both SI 
and SD interactions. As heavy nuclei are used in CDMS and XENON, the SD part in Eq. pS)) is several orders of 
magnitude smaller than the SI part, and thus the SD contribution may be omitted when comparing the results in 
those experiments. So, in this approximation, the WIMP-nucleus cross section for chiral interactions has the same 



form as that for vector interactions. The WIMP-nucleus cross sections in Eqs. (P5|) 
corresponding WIMP-nucleon cross sections as follows: 



Vector Int. 



Axialvector Int. 



Tensor Int. 



Chiral Int. 



'^V,xN 



<^A, xN 



<^T, xN 



^C, xN 





( Gv,N 








/ Ga,n 


tt{M^ + niNY 


\ V2 




f Gt,n 


Tr{M^ + TOAr)2 


[ ^/2 




/ Gv,N 




[ V2 



can be normalized to the 

(52) 
(53) 
(54) 
(55) 



Note that in Eq.(|55|) o'cxN is not the actual WIMP-nucleon cross section for chiral interactions, it just means a 
normalized quantity of the corresponding WIMP-nucleus cross section for a heavy nucleus A. In other words, a 
normalization procedure from Ea. (|i5)) to Eg. ([55]) is used: cfc^xN = [^^/(/^yi^^)] ' ^c,xA when A >• 1, with the 
reduced masses ijla = M^mAl (M^ + m^) and /ijv = M^mpf /{M^ + toat), where we have used the fact that Gv.p and 
Gv,n have the same order of magnitude: Gv,p ^ Gy.n- 

Eqs. (HH), ([5^ - ([551) predict the normalized WIMP-nucleon cross sections a^^^ for the different interactions as 
functions of the WIMP mass and the effective couplings to nucleons Gat, which further depend on the more 
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FIG. 3. The spin-independent (SI) WIMP-proton cross sections (dashed lines) for Dirac fermionic WIMPs with scalar, vector 
and chiral interactions. The results for chiral interactions are approximately identical to those for vector interactions. In the 
upper left and lower left frames, the results are shown for scalar and vector interactions for the case of universal couplings 
Gf = 10~^, 10~®, 10~^, 10~® and 10~^ GeV~^. In the remaining two frames, the results are shown for the case when the 
coupling constants proportion to the fermion mass m/, Gf x (1 GeV/m/) = 10~*, lO"'^, 10~^, 10~^ and 10~* GeV~^. Also 
shown as solid and dotted curves are the current upper bounds from the experiments of CDMS II (2009) [3l and XENON 100 
(2010) m. 



fundamental efFective couplings to quarks Gq. To illustrate the results and compare "with the experimental bounds, 
we draw the curves of a^^N "^s. for fixed coupling constants in these efFective models. 

In Fig.|3l we show the SI cross sections for a Dirac WIMP elastically scattering with a proton in the efFective models 
of scalar, vector and chiral interactions, respectively. The results for chiral interactions are approximately identical 
to those For vector interactions due to Eqs. ([5^ and ([55)) . so they are shown in the same Frame. We still consider 
the two types oF coupling constants. In the upper leFt and lower leFt Frames oF Fig. [31 the results are shown For scalar 
and vector interactions For the case oF universal couplings Gf = 10^^, 10^^, 10^^, 10~^ and 10~^ GeV~^. In the 
remaining two Frames, the results are shown For the case when the coupling constants proportional to the Fermion mass 
m/, Gf X (1 GeV/m/) = 10"^, 10""^, 10"^, 10"^ and 10"'' GeV"^. We also show in Fig.[2]the current upper limits 
For the SI WIMP-nucleon elastic scattering From the experiments oF CDMS II (2009) H and XENONIOO (2010) HH 
For comparison. 

In Fig. 01 wc show the SD cross sections For a Dirac WIMP elastiscally scattering with a neutron in the effective 
models oF axialvcctor and tensor interactions, respectively. In the upper leFt and lower leFt Frames oF Fig.jJl the results 
are shown For axialvector and tensor interactions For the case oF universal couplings Gf = 10~^, 10"'^, 10~^, 10~^ and 
10-2 GeV^^ In the remaining two Frames, the results are shown For the case when the coupling constants proportion 
to the Fermion mass m/, Gf x (1 GeV/m/) = 10"^, 10"'^, 10"^, 10~^ and 1 GeV"^. Also shown in Fig. S] are the 
upper bounds For the SD WIMP-neutron elastic scattering From the experiments oF CDMS (2005) [i^ and XENONIO 
(2008) [13. 
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FIG. 4. The spin-dependent (SD) "WIMP-neutron cross sections (dashed hnes) for Dirac fermionic WIMPs with axialvector 
and tensor interactions. In the upper left and lower left frames, the results are shown for axialvector and tensor interactions 
for the case of universal couplings Gf = 10"^ 10"^ lO""*, lO"'' and 10"^ GeV^^ In the remaining two frames, the results are 
shown for the case when the coupling constants proportion to the fermion mass m,/, Gf x (1 GeV/m/) = 10~*, 10~^, 10"'^, 
10^^ and 1 Ge'V"^. Also shown as solid and dotted curves are the upper bounds from the experiments of CDMS (2005) [i^ 
andXENONlO (2008) 



From Figs. [3] and [4l we see that the experimental constraints for the SI interactions (i.e., scalar or vector interactions) 
are much more stringent than those for the SD interactions (i.e., axialvector or tensor interactions). In addition, the 
experimental constraints in the case of universal couplings are more stringent than those in the case of Yuka^wa-like 
couplings {Gf (X mf) for both the SI and SD interactions. Due to the factor 4 difference in the WIMP-nucleon cross 
sections between us and that in [s^ our predicted WIMP-proton SI cross sections (dashed lines) in Figs. [3] arc lower 
than them. 

We summarize that the scalar, vector and chiral interactions are mainly constrained by the experimental SI upper 
limits of WIMP-nucleus elastic scattering, while the axialvector and tensor interactions are constrained by the SD 
upper limits of WIMP-nucleus elastic scattering. However, the other 6 types of interactions, £p, £sp, ^ps, ^va, ^av 
and Cf, cannot be constrained by direct detection experiments. 



V. INDIRECT DETECTION 



In addition to the direct search method for the WIMP DM at underground laboratories, an indirect detection 
method is used to look for the DM annihilation or decay products which include neutrinos, gamma rays, positrons 
and antiprotons. These particles can be detected by cosmic ray experiments. For the charged particles they are 
deflected by the Galactic magnetic field and interact with the interstellar medium when they propagate in the Galaxy. 
Therefore we have to study the propagation process of the charged particles to compare predictions with observations. 



13 



The propagation of cosmic rays in the Galaxy can be described by [56 



d_ 

dp 



Pi' - f (V • V,)i\ - —^p^-^ (56) 

3 i Tf Tr 



where ip = ip{r,p,t) is the number density of cosmic-ray particles per unit momentum interval, Q{r,p) is the source 
term, D^x is the spatial diffusion coefficient, Vc is the convection velocity, Dpp is the diffusion coefficient in momentum 
space describing the reacceleration process, p = dp/dt is the momentum loss rate, Tf and are the time scales for 
fragmentation and the radioactive decay, respectively. The most accurate method to treat the propagation is to solve 
Eq. ([55]) by a numerical code, GALPROP [S^H^l- For given source functions and boundary conditions, GALPROP 
can solve the equations for various cosmic-ray species and give reasonable fit to many cosmic ray data. 

Antiprotons are rarely produced in usual astrophysical processes, and the observed antiproton-to-proton flux ratio 
p/p is smaU, typically 0(10"^) ~ 0(10"^), from 100 McV to 100 GeV in kinetic energy. The annihilation of WIMPs, 
however, produces protons and antiprotons in equal numbers. Thus the p/p spectrum may be sensitive to the 
annihilation of WIMPs. In this work we use the antiproton-to-proton flux ratio p/p measured by the satellite-borne 
experiment PAMELA [s^] to constrain the effective models. 

The source term of antiprotons contributed by the annihilation of Dirac fcrmionic WIMPs is given by 



Qann (x, E) 



( dNp 
\dEp 



(57) 



with 



(o'annW)g, and Bq 



(o-ann'y)t 



-, q = u,d,s,c, b, t, 



(58) 



where p(x) is the DM mass density distribution of the Galaxy, {dNp/ dEp)^ is the number per unit energy interval 
of the antiprotons produced by the annihilation of a pair of WIMPs in the qq channel, and Bq is the branching ratio 
of the qq channel. The source term of protons is similar to Eg. ([57)) . We use the Monte Carlo program PYTHIA [H^ 
to simulate the particle production of the WIMP annihilation processes, and pick out the events in which the final 
states are p to build the {dNp/ dEp) ^ spectrum in each quark-antiquark channel. Due to the different quark masses, 
the {dNp/ dEp) ^ spectra in the six qq channels arc slightly different. 

The NEW profile [G^I is taken to describe the DM mass density distribution of the Galaxy: 



p{r) 



{r/rs) (1 -I- r/vs) 



(59) 



where Ps is the characteristic density, and is the scale radius. We choose Ps = 0.334 GeV/cm'^ and = 20 kpc. 
These values guarantee the local DM density p (8.33 kpc) = 0.4 GeV/cm^, which is consistent with the recent results, 
such as [m and [13 . The CDM particles in the Galaxy should obey the Maxwell-Boltzmann velocity distribution 
/ (fo) = (M^/27rfcBr)^^^ exp (— M^WQ/2fcsT) . Their velocity dispersion v = (vq) is chosen to be the canonical 
value 270 km/s d. Since the DM particles in the Galaxy nowadays are extremely nonrelativistic, i'm01 almost equals 
Wrei by neglecting the 0{v'^ / c^) term, and we need not to distinguish them. And Wroi has the same value in different 
frames with the extremely nonrelativistic limit. Thus, as a very good approximation, we expand (cTannf) to be 
(cannw) — a + blv'^) with the coefficients a and b obtained in Eqs. (P7)) - ([57)) . The relation between (i;^) and v is given 
by = (vf - 2vi • + vi) = 2 {vl) = 2v\ 

In the calculation of the p/p spectrum with GALPROP, we adopt the Galaxy propagation model with diffusion 
and convection, and set the half-height of the Galaxy propagation halo to be z^, = 4 kpc. Although the GALPROP 
expected p/p spectrum without the DM contribution can fit the PAMELA result well, as shown in the left frame 
of Fig. O the present data cannot rule out the possibility that the p/p spectrum might receive a small portion of 
contribution from the DM annihilation. If we add the DM contribution to the p/p spectrum, however, the total 
spectrum (GALPROP background -I- DM contribution) may deviate from the PAMELA result, as shown in the right 
frame of Fig.[5j For the universal couplings and the Yukawa-like couplings (G/ oc m/) in our numerical calculation, the 
DM contribution is actually determined by two parameters, ((Tannw) 
is, the smaller the DM contribution will be. We must let (era 
within an acceptable deviation range of the PAMELA result when we take into account the DM contribution. 

Now we calculate the value to set upper bound on the DM coupling constants from the p/p data. Since {(Ja,raiv)tot 
monotonously depends on the coupling constants G/, we can derive the 3cr upper bounds on G/ for fixed in the 



aim ^ 

/jQj and M^. For fixed M^, the smaller ((TannW/tot 
J be small enough to keep the total p/p spectrum 
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FIG. 5. The p/p spectrum calculated by GALPROP and that measured by PAMELA [5^ (error bars). In the left frame, we 
show the GALPROP predicted p/p spectrum without the DM contribution, where the Solar modulation potential $ is taken 
to be 335 MV to fit the PAMELA data. In the right frame, we show the predicted total p/p spectrum (GALPROP background 
+ DM contribution) in the effective model of scalar interactions as an example. = 1 TeV and ((Tannu)^^^ = 5.0 x 10^'* cm? /s 
are set here. 



cases of universal couplings and Yukawa-like couplings {Gf oc mf) for each type of effective interactions, as shown 
in Fig. [6l Here the solar modulation potential is set to be $ = 335 MV, which gives minimal for the background 
only. 

We see that several pairs of upper bound curves in Fig. [5] are nearly identical, which is similar to the situation of 
Fig. [T] The reason is the same as explain before. As pointed out in Sec. IIIH the quantity fiannW in each pair differ 
only by terms of O(w^) and/or terms of m'j/Al'^. Using the abbreviated notation defined in Sec. Illli we may call the 
pairs of nearly identical upper bound curves by S ~ SP, P ~ PS, V ~ VA, and T ~ C. These 4 pairs are the same 
as those in Fig. [T] in Sec. IIIIl The exception in this situation is that the upper bound curves for axialvector (A) and 
axialvector-vector (AV) interactions are rather different. It may comes from the fact that aAv. annV is of order 0{v'^), 
while cr A, annV is of Order 0(u°) and proportional to m'j. The difference of the curves for A and AV is somewhat 
like the difference of the cases of universal couplings and Yukawa-like couplings (Gf cx mj) for the same effective 
interaction. In addition, we note that in Fig. [5] the upper bound curves for scalar (S), scalar-pscudoscalar (SP) and 
axialvector-vector (AV) interactions lie well above the other curves, for their (JannV are of order 0{v'^). And we note 
that there are downward bends in the upper bound curves at about ~ nit ~ 171.2 GeV account for the X X ^ 1 1 
threshold effect. 



VI. VALIDITY REGION OF EFFECTIVE MODELS AND COMBINED CONSTRAINTS 

In this section, let us discuss the validity region where the method of effective theory can be used. For a generic 
4-fermion interaction operator x^ixf^^f , the mass dimension of the coupling Gf is —2. Since we have used 

the two types of coupling constants, i.e., the universal couplings and the Yukawa-like couplings {Gf oc mf), in the 
numerical calculation throughout the last three sections, let us consider them case by case: 

• For the universal couplings, we can write the coupling as ^ = where A is the cutoff energy scale and a 
is the coupling of the fundamental theory beyond A, which may be of order 1. The transfer momentum of the 
annihilation process X X ~^ f f niust be well below the cutoff, that is, 2M^ ^ A, so that the effective theory can 
be used. On the other hand, a weakly coupled UV completion of the effective theory usually requires a < Att 
such that the perturbation can be adopted [IiIIg^. From the above 3 relations, we obtain 

(60) 

• For the Yukawa-like couplings {Gf cx rrif), we have ^ = , where A is a cutoff and a is of order 1. Likewise, 
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FIG. 6. The 3(t upper bounds on the couphng constants Gf from the PAMELA p/p spectrum [5g] in each effective model 
of 4-fermion interaction operators. In the upper frame, resuhs are given for the case when the effective couphngs to aU the 
standard model fermions are equal (universal couplings). In the lower frame, results are shown for the case when the coupling 
constants are proportional to the fermion mass m/. 
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we still have 2My ^ A and a < An. From the above 3 relations, we obtain 



^ « ^ ■ (61) 



Eqs. (|60| and (|6T|) can be used to set the valid regions of the effective theories. Considering these validity conditions 
altogether with the other phcnomcnological constraints, we get the combined constraints on the effective models. In 
Figs.[7]-[T4l we show the combined constraints on the coupling constants G/ of Dirac fermionic WIMPs with scalar (S), 
pseudoscalar (P), vector (V), axialvector (A), tensor (T), axialvector-vector (AV), alternative tensor (T), and chiral 
(C) interactions, respectively. The constraints for the pseudoscalar-scalar (PS) interactions are nearly the same as 
those for the P interactions. The constraints for the scalar-pseudoscalar (SP) and vector-axialvector (VA) interactions 
are nearly the same as those for the S and V interactions, respectively, except for that SP and VA interactions are 
insensitive to direct detection experiments. Note that only the effective models of the S, V, A, T, and C interactions 
suffer constraints from direct detection experiments. The SI constraints on the S, V and C interactions arc much 
more stringent than the SD constraints on the A and T interactions. 

From these figures, we can get interesting results. For the scalar and vector interactions that induce SI scattering 
with nuclei, the constraints from direct detection can be much stronger than the other constraints for the universal 
couplings, but weaker if the coupling is Yukawa-like. This is easily understood by considering the relation of DM- 
nuclcon coupling constants and the DM-quark coupling constants given in Sec. IV. But for the axialvector and tensor 
operators that induce the SD interaction with nuclei we note that the direct detection constraints arc much weaker 
than the indirect p/p and relic density constraints. Further, we also note that in several cases the p/p constraints can 
be stronger than the relic density constraints for DM mass lighter than 100 GeV. If the direct or indirect detection 
constraints are stronger than relic density, the constraints on the effective coupling are so weak that the thermal 
production may overdose the universe. Therefore the DM models in such cases should be excluded, or else some 
exotic entropy generation process should occur after DM freeze out. 




FIG. 7. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with scalar (S) interactions from relic 
density, direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. The 
yellow region denotes invalid parameter space of effective field theory. The left frame is shown for the case of universal couplings, 
while the right frame for the case of Yukawa-like couplings (G/ oc m/). The constraints for scalar-pseudoscalar (SP) interactions 
from relic density, PAMELA p/p ratio, and validity of effective theory are nearly the same as above, but SP interactions are 
insensitive to direct detection experiments. 



VII. CONCLUSION 



In this work we give a general analysis of the 4-fermion interaction between the DM and the standard model 
particles. We have considered the most general form of the 4-fermion operators and corrected some errors in the 
previous works. We find that the constraints from DM relic density, DM direct detection and indirect detection oip/p 
data are complementary to each other. Generally, the SI constraints are the most stringent while the SD constraints 
are quite weak. For light DM (< 70 GcV) the p/p data give very strong constraints on the interaction. 
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FIG. 8. Combined constraints on coupling constants Gf of Dirac fermionic WIMPs with pseudoscalar (P) interactions from 
relic density, PAMELA p/p ratio, and validity of effective theory. The constraints for pseudoscalar-scalar (PS) interactions 
from relic density, PAMELA p/p ratio, and validity of effective theory are nearly the same as above. 




FIG. 9. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with vector (V) interactions from relic 
density, direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. The 
constraints for vector-axialvector (VA) interactions from relic density, PAMELA p/p ratio, and validity of effective theory are 
nearly the same as above, but VA interactions are insensitive to direct detection experiments. 



Assuming that one operator dominates the effective interaction between DM and the SM fermions, we find that 
some cases get so strong constraints that the universe will be overdosed by DM thermal production. In such cases 
the DM models are actually excluded assuming a standard cosmology. As a summary, in Tab. lU we indicate the 
excluded regions of given by direct and indirect experiments for Dirac fermionic WIMPs with various effective 
interactions. We find that recent direct detection experiments only exclude some regions of for the scalar, vector 
and chiral interactions with universal couplings, and for the scalar interaction with Gf oc mj. The PAMELA p/p 
spectrum, however, excludes some small Af^ regions (< 70 GeV) for most of the effective interactions. 
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FIG. 10. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with axialvector (A) interactions from 
relic density, direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of effective theory. 




FIG. 11. Combined constraints on coupling constants Gf of Dirac fermionic WIMPs with tensor (T) interactions from relic 
density, direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of effective theory. 




FIG. 12. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with axialvector- vector (AV) interactions 
from relic density, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 13. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with alternative tensor (T) interactions 
from relic density, PAMELA p/p ratio, and validity of effective theory. 




FIG. 14. Combined constraints on coupling constants G/ of Dirac fermionic WIMPs with chiral (C) interactions from relic 
density, direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. 
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